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We present a formalism for calculating the absorption coefficient of a pair of coaxial tubules. A 
spatially nonlocal, dynamical self-consistent field theory is obtained by calculating the electrostatic 
potential produced by the charge density fluctuations as well as the external electric field. There are 
peaks in the absorption spectrum arising from plasma excitations corresponding either to plasmon 

t-H or particle-hole modes, fn this paper, we numerically calculate the plasmon contribution to the 

absorption spectrum when an external electric fiefd is applied. The number of peaks depends on the 
radius of the inner as well as outer tubule. The height of each peak is determined by the plasmon 
wavelength and energy. For a chosen wave number, the most energetic plasmon has the highest peak 
corresponding to the largest oscillator strength of the excited modes. Some of the low-frequency 
plasmon modes have such weak coupling to an external electric field that they are not seen on the 

<^ same scale as the modes with larger energy of excitation. We plot the peak positions of the plasmon 

_i_ excitations for a pair of coaxial tubules. The coupled modes on the two tubules are split by the 

I Coulomb interaction. The energies of the two highest plasmon branches increase with the radius of 

the outer tubule. On the contrary, the lowest modes decrease in energy as this radius is increased. 

i i No effects due to inter-tubule hopping are included in these calculations. 

PACS numbers: 73.20.Mf, 71.45.Gm, 61.46.-w 
<D I. INTRODUCTION 

a 

^_i Recently, there has been a considerable amount of effort to understand the spectroscopic properties of carbon 

nanotubesP3E3 l n spectroscopic experiments, an external probe consisting of a beam of, for example, photons or 
electrons, perturbs the structure. Consequently, the ensuing charge density perturbation may be sustained by either 
single-particle excitations or collective plasmon modes. The density response may thus be measured through electron 
energy loss, photoemission, or absorption experiments, for example. So far, both theorists and experimentalists have 
investigated the plasma excitations of nanotubes using electron energy loss spectroscopy ( EELS) P"^ The theoretical 
formulation of EELS for nanotubes is now well established and several predictions concerning plasmon instabilities for 
nanotubes have been put forward. 9 In this paper, we present a theory for the absorption coefficient for a pair of coaxial 
tubules. The procedure can be generalized to multi-wall nanotubes. Several authoriP"^ have reported results for the 
optical absorption spectra of single- wall carbon nanotubes. The role played by the Coulomb interaction enters through 

r^. the dielectric function. This was obtained for single-wall nanotubes using the random-phase approximation (RPA) as 
well as time-dependent density functional theoryPEl \y e separate the contributions to the absorption coefficient due 
to plasmon and single-particle excitations by carefully examining the dielectric function. 

Our formalism relies on the use of the nonlocal dynamical response function. The reason is that this quantity 

l/~) determines the normal modes of oscillation which are identified as resonances. It is also responsible for the dielectric 
screening. We are able to separate the contributions from single-particle and plasmon excitations.^ For simplicity, 
we use a continuum model to formulate our theory for the absorption coefficient. However, the formula which we 
develop may be adapted to nanotubes with bandstructures obtained using a basis set of localized orbitalsP We employ 
^- our results for the single-particle eigenstates to determine the frequency-dependent nonlocal dielectric function. We 
neglect disorder by assuming that we are dealing with high-quality nanotubes. For single-wall carbon nanotubes, the 
mean-free paths typically exceed 1/j.m. We use Maxwell's equations for the electrostatic potential to calculate the 
dielectric response function of an electron gas on a nanotube in the RPAP Our analysis is specific to a double-wall 
nanotube but it can be extended in a straightforward way to a multi-wall nanotubesP 

Zhang, et alP^O calculated the dielectric response function for single- wall carbon nanotubes. These authors included 
the effect due to weak impurity scattering. Using their dielectric response function and our formula derived below, 
we would be able to identify the plasmon and single-particle excitations for multi-wall nanotubes with both electron- 
electron and electron-impurity scattering. Recently, the effects due to inter-tubule atomic hoppings and inter-tubule 
electron-electron interactions were included in calculations of the low-frequency electronic excitations of double-wall 
armchair carbon nanotubes!^ It was demonstrated in RefP^that the inter-tubule hopping significantly modifies the 
low-frequency plasmon excitations. For example, there are additional single-particle excitation channels and tunneling 
plasmon modes. Wc do include these effects arising from tunneling. However, within the framework of our calculations, 



these effects can be taken into account through the eiectron energy bandstructure. 

The rest of this paper is organized as follows. In Sec. [TTJ, we use a self-consistent field theory involving Poisson's 
equation and linear response of the charge density fluctuations to determine the induced potential. This potential 
is employed in the absorption coefficient whose resonances are determined by the plasma dispersion formula. We 
must separate the contributions from the single-particle and plasmon excitations when calculating the absorption 
coefficient. Gaps in the absorption coefficient for single-particle excitations correspond to pockets in the dispersion 
relation where the plasmons are not Landau damped. Numerical results demonstrating these assertions are presented 
in Sec. |III| Concluding remarks are given in Sec. |IV| 



II. GENERAL FORMULATION OF THE PROBLEM 

We now present a calculation for the absorption coefficient when an incident beam of light of frequency ui interacts 
weakly with a pair of coaxial tubules. In this case, linear response theory is applicable for calculating the rate at 
which energy is absorbed by the quantum system. The formalism we present below may be generalized to multi- 
wall cylindrical nanotubes consisting of an arbitrary number of coaxial tubules and to an array of parallel nanotubes, 
embedded in a material with dielectric constant €(,. We use Fermi's golden rule for the rate at which energy is absorbed 
and then apply the fluctuation-dissipation theorem to the density-density correlation function. 21 After some algebra, 
we find that the absorption coefficient is given in terms of the real (3?e) and imaginary (9m) parts of the Lorentz 
ratio a?Lorentz(w). This is defined as the ratio of the Fourier coefficient of the absorbed energy at frequency a; of a 
probing field to the square of its amplitude and we have 2 ^ 
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where p „„ (w) = [exp(fiu)/kBT) — 1] is the photon distribution function, the refractive index is 
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Here, <^ ex t ( r ) is an arbitrary external potential. However, for a uniform electric field, we have <^ ex t ( r ) = ~~ r ' E' 



ext 



Also, the electron density fluctuation 5n(r, uS) 



for an electric field in the direction with unit vector eo = E ext / E ext 

on the two tubules of inner radius R± and outer radius Ri arises from both the external electric field as well as the 
electric field from the induced potential, i.e., in linear response theory within the RPA, we have 
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Sn(r,u) =-el dr' ^ X ° (r,r';w) * tot (r',w) 



(4) 



where $tot ( r i w ) = ^ext ( r ) + ^ind ( r ' w ) 1S * ne sum °^ ^ ne external potential and the induced potential </> mc [ (r, uS) 
The response function for the j tubule is given by 
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In this notation, /o(e) is the Fermi-Dirac distribution function, ej tV is an energy eigenvalue on the j tubule and 
^ j v {y) the corresponding eigenfunction. If an electron with effective mass m* is confined on the surface of a tubule 
of radius Rj , whose axis is taken to be on the z axis, the eigenf unctions and eigenenergies are 



*iAw) = ^p^^M^r^ >*;M = s (r Rj) , (6) 
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with v — {k z ,l}, k z = — — , n,l = 0, ±1,±2, ••• and L z is a normalization length. Plasmon excitations can be 

z 

dp 
obtained by calculating the density matrix from its equation of motion ifi-pr- = [H,p\. For small perturbations, 

p = Po + Pi and H = Hq — e</>, where po is the equilibrium density matrix and p± its perturbation, Hq is the 
unperturbed Hamiltonian and (j>(r, ip, z; t) is the fluctuation in the electric potential corresponding to p\. 
After we substitute Eq. ^ into Eq. Q, we obtain the following result for Eq. ^, i.e., 

CkLorontz(w) = -e 2 ^ ^j,vv>{ U ) {Qj,uu> +Uj, vv >{u))} Mjy u , (8) 

where 

J |V</>ext( r )l 



Q ilVV , = J dv ext (r) **»%„,(!•) (10) 

are transition matrix elements and we have also introduced 



CWM = J dr ind (r,w) *^(r)* Av /(r) . (11) 

In the presence of an external potential </> ex t > the induced potential is a solution of Poisson's equation 



V 2 ind (r,w) = -^ fe M (r,w)-e / dr '5I X°(r,r';a;) ^ext( r ') 



(12) 



The electrostatic potential <fo n d( r ' ^> z > w ) can ^ e Fourier transformed with respect to the variables ip and z to give 
md (r, p, z; w) = X] fc.(r, g„ w)e^ V 1 "" , q, = ^ , n, L = 0, ±1, ±2, • • • , (13) 
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so that Poisson's equation becomes 



3=1 

Here. 
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and pi denotes the perturbed density matrix to lowest order in the external perturbation. The polarization function 
of the electron gas on each nanotube is defined by 



Xj,l(<7 z ,w) 



l.l' = -c 



2 E 

1,1'-- 
OO 

E 



dk z /o(ej,fc,,;) - fo{tj,k z -q z ,i-L,) 

2-K fiUJ + €j,k z -q z ,l-L - Cj,k M ,l 



dk 



tU) 



2 rrO 



m(k z ,l;k z ~q z ,l-L) 



(16) 



Also, in Eq. (14), we have introduced the quantity B^'(q z ,u>) which depends on the external electric field as follows 
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We now write the solution of Poisson's equation for 0i(r, q 2 ,o;) in Eq. (14) as 



f_ Ci/ L (^r), r<#i 

^(r,fe,w) = <^ CalUqzr) + C 4 K L (q z r), Ri < r < R 2 , 
I C 2 K L (q z r), v > i? 2 



(18) 



where II (x) and if/, (x) are Bessel functions of imaginary argument and C±, C 2 , C3 and C4 are functions of q z and w 
but are independent of the radial coordinate r. These coefficients are determined from the co ntinuity of 4>l{t, q z ',uj) 



at r = i?i,i?2 and the step- like change of — — at r = Ri,R 2 , obtained by integrating Eq. (14) across the surfaces of 

or 
the tubules. These boundary conditions together yield 
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In this notation, the prime on Il(%) and Kl(x) indicates that a derivative with respect to the argument must be 
taken. This set of simultaneous equations show that C±, C 2l C3 and C4 are determined by the induced potential 
through A j'{q z ,uj), B^(q z ,uj) on the right-hand side. On the other hand, the induced potential is determined by 
Cj in Eq. (18). This self-consistent field theory allows us to obtain the following analytic results for the coefficients 
as 
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.(i)/ 



i<n< 



X A^{q z ,w) + B^{q z ,uj) -R 2 K L {q z R 2 ) A? (&,«) + B^'{g s ,u) 



I W/ 



>( 2 )/ 



(20) 



C 2 (q z ,u) = -R 1 I L (q z R x ) 
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C 3 (q z ,uj) = -R2K L {q z R 2 ) A^>(q z ,oj) + B^(q z ,uj) 



|( 2 ), 
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(22) 



C i {q z ,u) = -R l I L (q z R 1 ) A£>(q z , to) + B ( l\q z , u) 



(23) 



Making use of these results in Eq. (18), we obtain c/) L (Ri,q z ,uj) — CilL(q z Ri) and 4>L(R 2 ,q z ,uj) — C 2 K L (q z R 2 ) on 
the surfaces of the tubules. Combining these results for <f>i,(Rj,q z ,oj) on the surface of the tubule with Eq. ( 15 1 where 

Af* '(q z ,uj) is given in terms of </> L (-Rj, q z , w), these self-consistent field results yield 
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V L (Ri,R 2 ;q z ,uj) = e\ j '{q z ,u))e\ j '(q z ,u) - V 12 (q z )xi,h(q z , w)x2x(g 2 , w) 



In the absence of an external perturbation, B^'(q z ,uj) = and the nontrivial solutions for .A^' (g z ,w) in Eqs. (24) 



(j). 



(26) 



and (25) correspond to V^(Ri, R 2 ; q Zl w) = 

We may now express M^ vv i in Eq. ([9]) and Uj mv i{lo) in Eq. (11) in terms of Fourier series expansions with respect 



to the variables <p and z. This was done for the induced potential in Eq. (13). We then have 
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where X (r) = 4>ey± ( r ) / l^^ext ( r )l ■ We a ^ so use similar Fourier expansions for Qj >vv ,' and Uj,w' in Eqs. (10) and 
(11), respectively. When the eigenfunctions are substituted into Eq. (27), then we find that the only non-zero matrix 
elements occur when V = I + L and k' z = k z + q z . 

We can separate the contributions to the absorption coefficient from plasmons and single-particle excitations. The 
normal modes of oscillation for each tubule when the coupling V\ 2 between them is neglected are obtained by solving 
££* (q z , w) = 0. Below, we derive explicit results for the induced potential for a single tubule of radius R. 

For a single tubule of radius R, we write the Fourier component of the induced potential as 



(fa(r,q z ,u) 



_ J C x I L {q z r), r < R 
C 2 K L (q z r), r>R 



(28) 



After applying the boundary conditions on the surface of the tubule to the continuity of the potential and the 
discontinuity of its derivative, as seen in (14), we obtain 



Ci(q x ,u) = -[A L (q z ,oj) + B h {q z ,Lj)}RK L (q z R) 

C 2 (q z ,io) = -[A^q z ,u)+B h (q z ,u)]RI L (q z R) , (29) 



where A^ and B^ are as denned above for a tubule of radius R. We now use the results in Eq. ( 29 ) to obtain the 
electrostatic potential 4>l{R, Izi^) on the surface of the tubule. When the derived result is substituted into Eq. ( 15 1, 
we obtain 

A L (q z ,u;) + B L (q z ,u J )= BL l qz > 0j) , (30) 



which can now be employed to calculate <j>i,(r, q z ,w) in all regions by making use of Eqs. ( 28 ) and ( 29 ). Here, e-^(q zi w) 



is given by s^'(q z ,uj)Eq. (26 1 for a tubule of radius Rj = R. 



In the absence of an external electric field, we must set </> e xt( r ) m Eq- ( 17) equal to zero. This results in B^ = so 
that Eq. (30) reduces to ei J (q z ,uj)Ai J (q z ,uj) — 0, where ei,{q z ,u)) — l + {e 2 /e^)lL(q z R)Ki J {q z R)xi J {qz,u) is the dielectric 



function for the nanotube. A nontrivial solution for Al(<7 z , u>) leads to the plasma dispersion equation ei,(q z ,uj) = for 



the self-sustaining plasma oscillations. However, in the presence of an external electric field, the solution of Eq. ( 30 ) for 



Ai,(q z ,uj) in terms of B^(q z ,uj) leads to the coefficients C\ and C 2 in Eq. (29) as Ci = —RKp{q z R)B^{q z ,u)/t\J^ z ,u}) 



and C 2 = —RlL(qzR)Bh(q z ,u)/eL(qz,w). When these results are substituted into Eq. ([28]), we obtain a closed 



form analytic result for 4>jj(r, q z ,uj) and consequently ^>j n( j(r, (p, z; u) in Eq. (13) by making use of the Fourier series 
expansion. We may now use these results to obtain U aa i (u) in Eq. ( |11| which could be employed in the Lorentz ratio 
in Eq. Q for calculating the absorption coefficient /3 a b s (w) in Eq. (TTF! 

We now turn to analyzing the contributions to aLorcntz(w) in Eq. ([8]). This means that we must examine when 
the factor 9m x?r( a; ) ( / , L(-Ri,<3'z7 w ) ¥" 0- Thus, we must determine 9m T>i J (Ri 7 R 2 ;q z ,uj)~ 1 = - „ 2 ' 2 where 2? K 
and T>\ are the real and imaginary parts of 2?l- Therefore, the only non-zero contributions to /3 a h s (w) arise when 
either T>\ is finite, i.e., from the single-particle excitations, or when both 2?r and T>\ vanish simultaneously. This latter 
contribution comes from the plasmon excitations which are not Landau damped by the particle-hole modes. 

III. NUMERICAL RESULTS 

In Fig. [2] we plot the plasmon contribution to the absorption coefficient when L = and q z = 0.15 kp as a function 
of the incident photon energy for inner radius R\ = 10.0 A and several values of the radius R 2 of the outer tubule 
for a pair of coaxial tubules. In Fig. |3j the peak positions of the two highest branches of plasmon excitation energies 
for a pair of coaxial tubules are plotted as a function of R 2 /Ri for fixed i?i = 10.0 A. The parameters used to do 
our numerical calculations were chosen as follows. The electron effective mass m* — 0.25 m e , where m e is the free 
electron mass, the Fermi energy Ep = 0.6 eV and the background dielectric constant tb = 2.4POur calculations show 
that the number of peaks depends on the chosen wave vector q as well as the ratio R 2 /R±. However, for the range 
< q/kp < 0.2, the two highest branches are excited. There are some small amplitude oscillations on the curve in 
Fig. [3] representing the highest plasmon excitation energy. This behavior is caused by the variation in the number of 
occupied energy levels on the outer cylinder as R 2 is changed for a fixed number of electrons, i.e., Ep. The oscillations 
of the second highest branch are negligible because the pocket in the single-particle excitation spectrum where that 
branch exists is very narrow. The intensity (oscillator strength) of each peak is determined by the plasmon excitation 
energy. Consistently, the highest peak has the largest oscillator strength.^ Some of the less energetic plasmon modes 
have such weak coupling to an external electric field that they are not seen on the scale we used in Fig. [2J 

In Fig. |4j we plot the particle-hole mode contribution to the absorption for a pair of coaxial tubules. The curve in 
Fig. HI is discontinuous Since the imaginary part of T> is zero in the pockets where the plasmon frequencies exist. The 
gaps on the w-axis correspond to the pockets in the excitation spectrum where the collective plasmon modes exist. 
The dispersion of plasmon and single-particle excitations is given in Ref. — citenew6. , this accounts for the gaps on 
the frequency axis of Fig. [4] We chose q z — 0.15 kp and all other parameters for the electron effective mass, Fermi 
energy and background dielectric constant are the same as Fig. [2] Unlike the plasmon contribution in Fig. [2j the 
particle-hole modes contribute over a wider range of frequency. The plasmon excitations occur within pockets between 
the particle-hole modes where they are Landau dampedP In Fig. pi we plot these particle-hole mode frequencies 
versus q z /kp for the pair of coaxial tubules used in Fig. HI Numerical results for plasma excitations between subbands 
for which L = 1, 2, • • • may also be obtained from our formula for the absorption coefficient. 



IV. CONCLUDING REMARKS AND SUMMARY 

In this paper, we gave a comprehensive formalism for calculating the absorption coefficient for a double-wall nan- 
otube. The calculation was carried out using a self-consistent field approach for the induced potential and charge 
density fluctuations on the tubules. Our result is given in terms of the electron-electron interaction on each tubule 
and between the two tubules. We did not include any effects due to inter-tubule atomic hoppingPS The effective 
dielectric function 2? L (Ri , R2 ', Qz , oj) for the pair of tubules is expressed in terms of the dielectric functions e^ (q z , oj) 



and £l {q z ,uj) for each tubule, as shown in Eq. (26 1. The electron gas model for each tubule was used to simplify the 
calculations. However, a more realistic model whose energy bands are obtained using a tight-binding approximation 
for example, would be incorporated through the susceptibility Xi,l(^,w) and X2,l(<7z,w) for each tubule. We showed 
that the loss function Im [l/T>^(Ri, R 2 ; q z ,ui)] can be separated into contributions due to plasmon excitations and 
particle-hole modes. Figures [4] and [5] show that there are pockets within the particle-hole continuum where there 
is no Landau damping of the collective plasmon excitations. These regions could only be determined by separating 
the contributions to the loss function from the plasmon and single-particle excitations, as we described above. This 
separation would allow direct comparison between theory and experimental results of the absorption spectrum for 
plasma excitations on nanotubes. 
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FIG. 1: Schematic of a pair of coaxial tubules. 
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FIG. 2: The L — plasmon contribution to the absorption coefficient (arbitrary units) versus photon energy for q z = 0.15 kp, 
where Uf — y2m*£jf//! 2 with the radius of the inner tubule R\ — 10.0 A and several values of R2. For clarity, the curves are 
shifted. From the bottom, the curves correspond to R2 = 20 A, 30 A and 40 A. The parameters used in the calculation are 
given in the text for m* , Ef and e^. 
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FIG. 3: Peak positions for the L — and q z — 0.15 Uf plasmon excitations for the absorption coefficient for a pair of coaxial 
tubules as a function of R2/R1, for inner radius Ri — 10.0 A and outer radius ife- The parameters used in the calculation are 
given in the text for ra* , Ef and tb and are the same as Fig. [2] 
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FIG. 4: The L = particle-hole mode contribution to the absorption coefficient (arbitrary units) versus photon energy for 

q z — 0.15 kjp, where Isf = \/2m*EF/Ti 2 ■ The radius of the inner tubule Ri = 10.0 A and R2 = 20.0 A is the radius of the outer 
tubule . Only the frequency range where there are particle-hole modes contributes. The gaps on the frequency axis correspond 
to pockets where plasmon excitations are not Landau damped. The parameters used in the calculation are given in the text 
for m*, Ef and et- 
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FIG. 5: The dispersion relation for the L — continuum of particle-hole modes (shaded regions) for the pair of coaxial tubules 
used in Fig. [4] 



